MEAN FIELD LIMIT FOR DISORDERED DIFFUSIONS WITH 
SINGULAR INTERACTIONS 
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Abstract. Motivated by consideration from neuroscience (macroscopic behavior of large 
ensembles of interacting neurons), we consider a population of mean field interacting dif- 
fusions in R m in the presence of a random environment and with spatial extension: each 
diffusion is attached to one site of the lattice Z d and the interaction between two dif- 
fusions is attenuated by a spatial weight that depends on their positions. For a general 
class of singular weights (including the case already considered in the physical litera- 
ture when interactions obey to a power law of parameter < a < d), we address the 
convergence as N — > oo of the empirical measure of the diffusions to the solution of a 
deterministic McKean-Vlasov equation and prove well-posedness of this equation, even 
in the degenerate case without noise. We provide also precise estimates of the speed of 
this convergence, in terms of an appropriate weighted Wasserstein distance, exhibiting 
in particular nontrivial fluctuations in the power law case when | ^ a < d. The example 
of FitzHugh-Nagumo oscillators is in particular considered. 
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1. Introduction 

The purpose of this paper is to provide a general convergence result for the empirical 
distribution of spatially extended networks of mean field coupled diffusions in a random 
environment. We are in particular interested in diffusions modeling the spiking activity of 
neurons in a noisy environment. To motivate the mathematical model we want to work 
with, let us consider, as a particular example, a family of stochastic FitzHugh-Nagumo 
neurons 

dVi (i) = (Vi (t) - Mf- - Wi (t) + j) dt + a v &Bj (t), 

dwi(t) = (aiibiViit) - Wi(t)fj dt + a w dBf{t), ' ' ' 

for % = —N,—N + 1, . . . , N — 1,N, with exterior input current I. The variable Vi(t) 
denotes the voltage activity of the neuron and Wi(t) plays the role of a recovery variable. 
[bJ {t) , Bf (t)) are independent Brownian motions modeling exterior stochastic forces. 
Depending on the parameters (aj,6j) e R 2 , the neurons exhibit an oscillatory, excitable or 
inhibitory behavior. Suppose that the precise values of u)i = (a,i,bi) are unknown, which 
will always be the case in a real-world application, but rather are given as independent and 
identically distributed random variables. From a point of view from statistical physics, 
this additional randomness in (jl.ip may be considered as a disorder. For simplicity we 
suppose that the uji are independent of the time t. 
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Equation can be written as 

ddi(t) = c(9 i ,uj i )dt + a- dBi(t), i> 0, i = -N, . . . , N, (1.2) 

using the shorthand notations 9 = (V, w), uj = (a, 6), c(6,ui) = (v — ^ — w + I,a(bV — 

w)j, B = [B v ,B w ) and a = y 7 ^ ^ V We suppose that the individual neurons are 

coupled with the help of a possibly nonlinear and random coupling term T (8i,Ui,dj,ujj), 
—N sj i,j sj N modeling electrical synapses between the neurons. The coupling intensity 
between neurons i and j will depend in addition on some weight (vE'jv(v) m &y 

be thought as a distance function, but not necessarily), so that the resulting system gets 
the following type: 



d0i(t) = c{6 i {t),oj l )dt 

N 

-X V. 

j=-N 

(1.3) 



1 ^ 



Particular examples of the weight are the P- nearest neighbor model ^jsf(i,j) = 1 if 
\i — j\ ^ P and otherwise, or a polynomial decay = \i — j\ a of the coupling 

strength. Note that higher-dimensional analogs of (]1.3p will also be considered in this 
paper. 

In the case of = 1> systems of the type (I1.3P are called mean field models (or 
weakly interacting diffusions) in statistical physics and have attracted much attention in 
the past years (see [24j ESI [32j HD]), since they are capable of modeling complex dy- 
namical behavior of various types of real-world models from physics to biology, like e.g. 
synchronization of large populations of individuals, collective behavior of social insects, 
emergence of synchrony in neural networks ([2j 134} lllj). and providing particle approx- 
imations for various nonlinear PDEs appearing in physics ([Zl HI [22], [5]). The most 
prominent example of such models is the Kuramoto model ([T| \21\ [3] 116]): 

K N 

d9i(t) = ujidt + — Yi sin (0, -0*) dt + adB^t), t ^ 0, i = -N,...N. (1.4) 

2 + 1 j=-N 

where K ^ is the intensity of interaction and 0j e S := ~R/2n. In the context of weighted 
interactions, a notable attempt to go beyond pure mean field interactions has been to 
consider moderately interacting diffusions (see [271 [23 EH] ) . 

All the statistical information of the neural ensemble is contained in its empirical dis- 
tribution 

1 N 

u { t N \de, du) : = ^ <W),«*)( d(? > duj ^ f > ( L5 ) 

j=-N 

that can be seen as a random probability measure. Since we are interested in the collective 
behavior of a large numbers of neurons, as it is the case for neural ensembles in the brain of 
living organisms, understanding the asymptotic behavior of ' as N — > oo is important. 

Under mild assumptions on ^jv, we will prove, as part of our main results in this paper, 
that v\ N ' converges to a deterministic measure ft(d0, duj) = qt(9,ui) d0/z(dw) where qt is 
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a weak solution of the McKean-Vlasov equation 

d t q t = ^dw e (aa T Veqt) - div e ^ jc(0,w) + J T(9,u>,9,Lj)q t (9,u>) d0/x(dw)j) . (1.6) 

A priori existence of weak solutions of (|1.6jl is obtained under mild assumptions in Propo- 
sition [22] below. Note that, of course, for degenerate noise, one cannot expect to have a 
density. The measure vt is called the mean field limit and we will provide in Theorem 12. 121 
(resp. Theorem 12. 1TH explicit rates on the convergence of v\ against vt in terms of an 
appropriate weighted Wasserstein distance. 

Note that the method followed in this paper significantly differs from the ones usually 
used in this context (i.e. propagation of chaos estimates, [211 [191 [26]). I* 1 that sense, this 
work should also be relevant even for the non disordered case without spatial extension. 



2. Mathematical set-up and main results 

2.1. The model. Fix JV ^ 1, T > and let be the hypercube [-JV, . . . , Nj d a Z d 
and |Ajv| = (2iV + l) d be its volume. We consider |Ajy| diffusions on [0,T] with values in 
the state spac^H X := R m for a certain m ^ 1. 

Each diffusion 9i is attached to the site i of A^r- The local dynamics of 9i is governed by 
the following stochastic differential equation which is perturbed by a random environment 
represented by a vector u>j e £ := R n (n ^ 1). 

d0i(t) = c(9 i ,u i )dt + a- dB^t), s^t s^T,ie A N . (2.1) 

where a e J{, mxm i s the covariance matrix, c(-, •) is a function from X x £ to X and (_Bj) 
is a given sequence of independent Brownian motions in X. The vectors {uJi)i<=\ N are 
supposed to be i.i.d. realizations of a law fj, and are hence seen as a random environment 
for the diffusions. 

When connected to the others, the diffusions interact in a mean field way with spatial 
extension: 

d9 i (t) = c(9i,ui i )dt+ T -— ; V r(9 i ,oj i ,9 j ,oo j )^ N (i,j)dt + a- dBi(t), 

0^t^T,ieA N , (2.2) 

where T is a function from (X x£) 2 to X, and for fixed N, is a function from Z d x Z d to 
[0, oo). The required assumptions for the function will be made precise in Section |2, 2. 41 
below. One should notice at this point that does not need to depend on the difference 
i-j- 

We suppose that, at time t = 0, the variables (6i(0))i ^ n are independent and 
identically distributed according to a probability distribution C(d9) on X. 

2.2. Examples of weight functions ^tv. The class of weights we address includes non- 
trivial space-dependences already studied in the physical literature. Interesting cases are: 



Note that it is also possible to choose X as the circle S in the case of the Kuramoto model, but we 
will stick to X := R m for simplicity. 
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2.2.1. The P-nearest neighbor model: this model, originally introduced and studied in 
dimension d = 1 (see [281 129] ) concerns the case where each diffusion 6i only interacts with 
a number P (smaller than N) of its nearest neighbors on Z. This model can be easily 
generalized to any arbitrary dimension d by 

d0i(t) =c(9i,u}i)dt+jj—r J] r(^,^,^, Wj ) df + a- dBi(t). (2.3) 

Here, the oscillator 6i e A^r only interacts with its neighbors within the box of size P. We 
address in this work the behavior of the system (|2 . 3[) in the case where P is proportional 
to N, that is 

P = RN, (2.4) 

for a fixed proportion R e (0, 1]. 

Remark 2.1. The case of R = 1 corresponds to the mean field case. Understanding the 
behavior of the system (|2.3f) in f/ie case of a pure local interaction (that is when P « N) 
does not enter into the scope of this work. In particular, we will not address the question 
of P of order smaller than N (e.g. P = RN a for some a < 1), whose behavior as N — > oo 
seems to be quite different. 

Under the assumption (I2.4p . the P- nearest-neighbor model (I2.3P enters into the frame- 
work of ()2.2p for the following choice of ^jy: 

Vi, j 6 A N , * N (i,j) := ^JL-l^, f^pj := XR (^pj . (2.5) 



2.2.2. The case of a power-law interaction: this model also considered in the physical 
litterature corresponds to the case where the interaction between two diffusions 9i and 9j 
decays as a function of || i — j \\~ a for some parameter a ^ 0. Note once again that the 
pure mean field case corresponds to a = 0. 

As observed in [9l [III [23l [31] on the basis of numerical simulations, it appears that 
the behavior of the system is strongly dependent on the value of the parameter a. The 
situation which is considered in this paper corresponds to the subcritical case where the 
parameter is smaller than the dimension: 



a < d. 



(2.6) 



The case of a ^ d is much more delicate and will be the object of a future work. We refer 
to § 12.2.31 below for further explanations on this case. 

As noticed in [T7], in the case of (|2.6p . one needs to renormalize the position i of the 
diffusion 6>j by a factor in order to have a proper continuous limit, that is 



*jv(i,j) := 



i 

2N 



3 

2N 



(2.7) 



Let us give an intuition of this renormalization: by standard arguments, it is easy to see 

j |p a is of order N d ~ a (since a < d). Conse- 



that the diverging series Yjjek 
quently, ttt-t^ 



2N 



IS of order J^jN 

• 00. 



nontrivial limit in (|2.2p . as 



d—c 



0(1), so that we should expect a 
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The renormalization by ^ maps Ajv = [—A", . . . , Nj d to a discrete subset of [ — i, |] 

converges to || x — y \\~ a which is an integrable function on [—5, ^] d - 
We refer to Section 12.2.41 for a precise definition of the continuous limit in this case. 

2.2.3. Some remarks on the supercritical case: the case of a ^ d is more delicate and 
requires more attention. Note that, to our knowledge, no proposition for any continuous 
limit has been made in the literature in this case. We are only aware of [9], where the 
system (j2.8|) is considered for finite N. 

One trivial observation is that the series YijeA N j& IP ~~ 3 IP" ^ s in * ms case already con- 
vergent. As a consequence, an interaction term of the form Yjjeh^ T Wj, 0j,0Jj) \\ i — j 

N j& 

in (|2.2p simply vanishes to as N — > 00. Hence, the correct model in this case is where 
the factor rJ— r is absent 

dB i (t)=c(6 i ,Ui)dt+ ^ i^^.p dt + ff' dBi(t). (2.8) 

The main difficulty for the derivation of the correct continuous limit in the case of (|2.8p 
lies in the fact that the interaction term XijeAjy F w i> w j) IP — j II ~° is n °t sufficiently 

mixing: if it exists, the McKean-Vlasov limit in this case should be random. We believe 
that the correct continuous limit should be governed by a stochastic partial differential 
equation instead of a deterministic PDE. This case is currently under investigation and 
will be the object of a future work. 

2.2.4. Periodic vs non-periodic boundary conditions. Instead of considering diffusions on 
Ajy, we can also suppose periodic boundary conditions, i.e. when An is replaced by 
Ajv,per := T^, where TV is the discrete A^-torus, that is [—A/", . . . , NJ with — A" and A" 
identified. The only thing that changes in what follows in the continuous model is that 
one should replace [ — i, |] by T d where T := [ — i, 2]/- I~A- Since the corresponding 
changes in the proofs of this paper remain marginal, we will restrict to the non periodic 
case and let the interested reader make the appropriate modifications in the periodic case. 

2.3. Notations and assumptions. 

Assumption 2.2 (Hypothesis on T and c). We make the following assumptions: 

• The function (9,lj) >—> c(6,cj) is supposed to be locally Lipschitz- continuous in 9 
(for fixed oj) and satisfy a one-sided Lipschitz condition w.r.t. the two variables 
(9,lj): 

V(p,<j),(B,u), (9-9,c{9,uj)-c{9,Cd)) ^L(J|0-£|| 2 + ||cj-w|| 2 ) , (2.9) 

for some constant L (not necessarily positive). We suppose also some polynomial 
bound about the function c: 

V(0,£j),||c(0,£j)|| < III c |||(1 + 11011" + (2.10) 

for some constant \\\ c\\\ > and where k 2 and 1 1. 

• The interaction term T is supposed to be bounded by \\ T and globally Lipschitz- 
continuous on {X x £ ) 2 , with a Lipschitz-norm \\ T \\ Lip . 



so that 



2N 



J 

2N 
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We also assume that for fixed 6, to, to, the functions 6 >—> c(6,u) and 9 >—> T(0,u),6,uj) are 
twice differentiate with continuous derivatives. 

Remark 2.3. Assumption \2.2\ is in particular satisfied for the FitzHugh-Nagumo case. 
One technical difficulty is the dynamics is not globally Lispchitz continuous. This will 
entail some technical complications in the following. 

Note also that the constant \\\ c \\\ mentioned in (|2.9p does not take part in the estimates 
of Sections [^] to [6j It only enters into account in Section^ 

Assumption 2.4 (Assumptions on fx and (). We suppose that the initial distribution C 
of 6 satisfies the following moment condition: 

f \\9\\ K ({d6) < oo, (2.11) 
Jx 

and that the law of the disorder [i satisfies the moment condition: 

J || u fn{du)) < oo, (2.12) 

where t is given by ()2.10|) . 



Assumption 2.5 (Assumptions on the weight ^). In order to cover the case of both the 
P '-nearest model and the power-law interaction introduced in § \2.2[ we suppose that either 
Hypothesis \(H1 ) or Hypothesis (H2)\ is true: 

(HI) P -nearest-neighbor: 

Vi,j e A N , * N (i,j) := [-^f) > (2-13) 



where the function \I/ is equal to \R (")> introduced in (12.51) . 
(H2) Power-law: 

e A N , $ N (i,j) := * (-^, , (2.14) 

where ^ is supposed to be a nonnegative function on [ — ^, IJ dxd suc h that the 
following properties are satisfied: 

Xi(^) := sup \\x-a\\ a ^(x,a) <co, (2.15) 

Z 2 (*):= sup m^)-ny,x)\dx <oQ (2ig) 
-,»e[-H] - l|X - yll " Q 

|| x — a || 27 ^(x, a) — || y — a || 27 ^(y, a) 

Z 3 (*) := sup - < oo, (2.17) 



II 27— ot 

r 1 n d \\ x — V\ 



for some parameters a e [0, d) and 7 chosen to be 



7 e [a, f ) 2/ a e [0, 1 
7 = | otherwise. 



(2.18) 



Remark 2.6. iVoie £/ia£ we cow/d /lave chosen simply 7 = | in any case. -Bui £/ns would 



have led to worse convergence rates than the ones that we obtain below in Theorem 2.1' 
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Of course, the main prototype for Hypothesis (H2) is when \]/ (x,y) = \\x — y ||~ a , for 
a < d (recall (12.70 ). But, the assumptions made in (H2) cover a larger class of examples: 
the reader may think of the general case of ^(x, y) := ip(x,y) \\ x — y \\~ a , for a bounded 
Lipschitz-continuous function ip. Note also that the case of bounded Lispchitz interactions 
is also captured (take a = 0). 

2.4. The empirical measure. Similar to (jl.5p . except that we have to take into account 
the renormalized positions of the diffusions, let us consider for fixed horizon T and time 
t e [0, T], the empirical measure uj. N ^\ 

v t N \d6, duj, dx) := ^(t),^,^ d#, du, dx), (2.19) 

j 

as a probability measure on X x £ x [- |,|] Here 

\ j 6 An- (2.20) 



j r 1 1^ 



2iV E 



2' 2 



2.5. The McKean-Vlasov equation. The convergence of the empirical measure at t = 
is clear: since (#i(0),£jj)i <c j <j at are i.i.d. random variables sampled according to ( (x) 
the initial empirical measure i/g ^ converges, as N — > oo, to 

i/ ( d0, dw, dx) := C( d9)fj,( duj) dx. (2.21) 

In addition, it is easy to see that if t i— ► converges for all t, its limit t >— *■ v t should 
satisfies the following nonlinear McKean-Vlasov equation 

dt {vt , f) = (ut , ^div e {aa T V e f) + V e f • c(-, 

Vt, Vef-jr(;;0,Lj)*(;x)u t (d9, did, dx)^> , (2.22) 

where is the weight function introduced in § 12.31 and (6,oj,x) >-> f(6,uj,x) is any 

bounded test function that is of class C 2 w.r.t. 9 with bounded derivatives. Note that we 
use here the usual duality notation (u , /} = $ / du for the integral of a test function / 
against a measure v. 

Remark 2.7. An important remark about a priori properties of (|2.22|) is the following: 
taking a test function f in (|2.22p that does not depend on 9 implies 

<*b, /> = <"*,/>, Vie[0,T]. 

In particular, the marginal distribution of (uj,x) w.r.t. the measure v% is independent oft 
and equal to d//(x) dx. This implies that, for a kernel ^ in the class of singular weight we 
consider here, it is always integrable against vt, for all t, since the function y >—> \\ x — y 
is integrable w.r.t. to the Lebesgue measure on [ — i, 

Moreover, since the function c is supposed to have a polynomial growth (recall (|2.10p ). 
one has to justify in particular the term (yt , V#/ • c(-, •)) in (|2.22p (the others are easily 
integrable). Thus, one should look for solutions t *—> vt having finite moment: for all 
te [0,7 |. \ Xx£ \\6\\ R \\uj\\ L u t {d9, dw, dx) < oo. 

In particular, well-posedness in (I2.22p will be addressed within the class of all measure- 
valued processes satisfying the properties mentioned above. 
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Formally integrating by parts in equation (|2.22|) . assuming the existence of density 
z^t(d0, dcj, dx) = qt(0,uj,x)d9fi(duj)dx, qt satisfies 

Mt = 2 div 6» ((Jcr T V e qt) - div e (q t (6,uj,x)c(9,u))) 

/ r \ ( 2 - 23 ) 

— divg I q t (9, u, x) T(9, u, 9, u)^(x, x)q t (9, u, x) d9[i{ du) dx J , t > 0, 

In the case where a is non degenerate, one can make this integration by parts rigorous: 
using the same arguments as in |15t Append. A], one can show that for any measure- valued 
initial condition in (|2.22p . by the regularizing properties of the heat kernel, the solution 
of (|2.22p has a regular density q t for all positive time that solves (|2.23|) . We refer to [T5J 
Prop. A.l] for further details. But of course, if a is degenerate, the strong formulation 
(|2.23p does not necessarily make sense and one has to restrict to the weak formulation 
(j!T2~2~j) in that case. 

2.6. Results. 

Proposition 2.8 (Existence of a solution to the McKean-Vlasov equation). Under the 
assumptions made in § \2.3l for any initial condition vo(d9, dw, dx) = £( d(9)//( du) dx, 
there exists a solution t>—*ut to f|2.22|) . 

The proof of Proposition 12.81 can be found in Section [3l Having proven the existence of 
at least one such solution in the general case, we specify now the problem to the case of 
Hypothesis |(H1)| (§ l2lTll and to the case of Hypothesis |(H2)| (§ 

2.6.1. The P -nearest-neighbor case: Suppose that the weight function *$>n satisfies Hy- 



pothesis (HI) 



Definition 2.9 (Test functions for P-nearest-neighbor). For fixed R e (0,1] and a e 
[ — 5> 2] 7 let Cr,o, be the set of functions f on X x £ x [ — i, i] of the form: 

f : (9, u), x) h-> g (8, u) -XR(x-a), 
where g is globally Lipschitz- continuous w.r.t. (8,uj): 

3C7 > 0,^(9, u, 9, cd), \g{e,u)-g{e,u)\ < C (\6 - 9\ + \u - ui\[) . (2.24) 

Let 

11 r 11 II 9{8,u) - g(0,Q) I! 

ll/lka : = SU P To — ITT^i =11 

e,e,w,ui || y _ y || + II ^ _ w II 

be the corresponding seminorm. 

Remark 2.10. Note that for any f e C Ra that is C l in the variable 9, the following 
estimate holds: 

W,u>,x, \\V e f(9,cj,x)\\ ^ \\f\\ R>a XR(x-a). (2.25) 
We now turn to the appropriate distance between two random measures: 
Definition 2.11 (Distance for P-nearest neighbor). For random probability measures A 
and v on X x £ x [ — i, i] , let 

1/2 



d R {\v) :=sup E|| </, A>-</, v)\ 
f v 

where the supremum is taken over all the functions f e Uoe[-l i] d ^R,a> su °h that \\ f \\ Ra ^ 1, 
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Our convergence result is given in the following 

Theorem 2.12 (Law of Large Numbers). Under Hypothesis \(H1 )\ for all R e (0, 1], we 

have: 

C 

sup d R {y[ N \v t ) < — - (2.26) 
where the constant C > only depends on T , Y , R and c. 

2.6.2. The case of the power-law interaction. Assume that the weight function ^> n satisfies 



Hypothesis (H2) 



In view of the form of \& in this case (recall Assumption 12 .5p . the main idea is to consider 
test functions (9,oo,x) t— * f{9,oj,x) that become regular when renormalized by || x — a \\ a . 
The seminorm || • || introduced in (|2.30p below should therefore be thought of as a weighted 
Holder seminorm. 

Definition 2.13 (Test functions for power-law interaction). For fixed a and 7 as in 
Assumption \2.5\ and for fixed a £ [ — \ , |] , let C a be the set of functions (9,u,x) <->■ 
f(9,to,x) on X x £ x [ — i, satisfying: 

• Regularity w.r.t. (9, 00): the function {9 , to) 1— >• \\x — a\\ a f(9,to,x) is globally Lipschitz- 
continuous on X x £, uniformly in x, that is 

3C> 0,^(9, to, 9, Co), \\x-af\\f(9,to,x)- f(8,td,x)\\ ^ C (\\9 - 9\\ + ||w - w||) (2.27) 

• Regularity w.r.t. x: the function \\ x — a \\ a f{9, to, x) is uniformly bounded 

1C > 0,\\x -a\\ a \\f(8,to,x)\\ ^C, (2.28) 

and the function x *— * \x — a| 27 f(9,uo,x) is globally (27 — a)-Hdlder, uniformly in 
(0,to): 

3C>0, Wx-af"' f{9,oo,x)-\\y-a\\ 2 " f f{8,oj,y) < C \\ x - y f^ a . (2.29) 
Denote by 

n,n \\x-a\\ a \f{8,uj,x)- f{8,uj,x)\ tt|l ,,. . „ 

ll/L := SU P 11. ^11 r =1 + sup || x - a (I \\ f(9,io,x)\\ 

9 - 9 + to - to e.ui.x 



\\ x - af 1 f{0,u,x) - \\y - a|| 27 f(9,io,y) 

+ SU P S i5F5 ( 2 - 3 °) 

$,u,x,y || x — y II ' 

the corresponding seminorm. 

Remark 2.14. Note that for any f e C a that is C 1 in the variable 9, the following holds: 

\/8,to,x, \\Vef(9,to,x)\\ < II/IIq ||Q . (2.31) 

|| x — a || 

The corresponding definition of the distance between two random measures is similar 
to Definition 12.111 given in the P-nearest neighbor case. The main difference here is that 
one needs to take care of test functions with singularities. Since those singularities happen 
at points of the form ^ (f° r some i and N) that are regularly distributed on [ — i , ^1 , 
we first need to introduce some further notations: for all integer K ^ 1, we denote by T>k 
the regular discretization of [ — \, |] with mesh of length 

d 





1 1- 


h 


-~ 2' 2. 



(2.32) 
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K > 1 



The appropriate distance between two random measures is then: 

Definition 2.15 (Distance for power-law interaction). Let a < d and p ^ 2 be defined by: 

V-{ 2 (2 33) 

where \x\ stands for the smallest integer strictly larger than x. On the set of random prob- 
ability measures on X x £ x [— i, let us define a sequence of distances \ d^(-, •) 
indexed by K ^ 1, between two elements A and v by 

4 ) (A,z,)=snp(E||</, A>-</,z,>n 1/p 
/ 

where the supremum is taken over all the functions f e (J aeV K , C a , such that \\ f \\ a ^ 1. 

i<: K' <: K 

Let us then define the distance d&, (■, •) 6y 

4 P) (A,^) := ^ 2 4 £: ^ (4 } (A^) a l) , (2.34) 

for a sufficiently large constant C (that depends on the parameters of our model). For a 
precise estimate on C , we refer to Proposition \6.5\ below. 

dp 

— CK~*}~ 

Apart from the weight e K2d — (which is precisely here to compensate the estimate 

that we find in Proposition 16.51 below), the definition of 0$ (•> ") exactly follows the usual 
Frechet construction (see e.g. |14|). 

Remark 2.16. The choice of the integer p in (|2.33p is made for integrability reasons that 
will become clear in the proof of Theorem \2.17\ One only has to notice here that p has 
been precisely defined so that its conjugate q (~ + ~ = 1) always satisfies qa < d. 

The main result of this work is the following: 

Theorem 2.17 (Law of Large Numbers in the power-law case). Under Hypo the sis \( H2 ) \ 

p. we have the following convergence: 



sup d%\v[ N \ 

^ t ^ T 





r i 








In AT 

d ) 






In AT 




I N d - a 



ifae [0,f), 

ifa = i (2.35) 
ifae (| ,d) , 



where the constant C > only depends on T, V, ^> , a and c. 

Note that the speed of convergence found in Theorem 12.171 is never smaller that iV~2 
which is the optimal speed for the case without spatial extension (recall the CLT results in 
the mean field case in [21]). Note also that, in the case where ^ a < rj, we have obtained 
a speed of convergence which is arbitrarily close to this optimal speed (since in that case 7 
is arbitrarily close to ^). We believe that the optimal speed in this case should be exactly 

, but the proof we propose in this work does not seem to reach this optimal result. 
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Nevertheless, in the case where we only consider a bounded Lispchitz- continuous weight 
function (i.e. with no singularity at all), the proof of Theorem 1 2 . 1 71 can be considerably 

simplified and one obtains a speed that is exactly N~z. 

Note also that the fluctuations when a e \^,d) appear to be nontrivial. A natural 
perspective of this work would be to prove a precise Central Limit Theorem in this case 
and to study the limiting fluctuation process in details. 

2.7. Well-posedness of the McKean-Vlasov equation. A straightforward corollary 
of Theorems [2T2] and ET7] is that uniqueness holds for the McKean-Vlasov equation (12. 22ft : 

Proposition 2.18 (Well-posedness of the McKean-Vlasov equation). Under the assump- 
tions made in § \2.3\ for every initial condition VQ(d6, da;, dx) = £(d#)//(du;) dx, there ex- 
ists a unique solution t >-* v t e M.\ ^C([0,T],<Y) x £ x [—5, 5] J to the McKean-Vlasov 
equation (|2.22[) . 



2.8. Organization of the remainder of the paper. The rest of the paper is organized 
as follows: Section [3] contains the proof of Proposition 12.81 concerning the existence of a 
solution to the McKean-Vlasov equation (|2.22p . Section [J] summarizes the main properties 
of the propagator of the nonlinear process associated with (I2.22|) . The proofs of the laws 
of large numbers are provided in Section for the P-nearest case (Theorem I2.12|) and in 
Section [6] for the power-law case (Theorem 12.17ft . An additional assumption of regularity 
is made from Section [H to El with is is discarded in Section [71 



3. The non-linear process and the existence of a continuous-limit 

The purpose of this paragraph is to prove Proposition 12.81 concerning the existence of a 
solution to the McKean-Vlasov equation fj2.22|) . This part is reminiscent of the techniques 
used by Sznitman (|33j) in order to prove propagation of chaos for non disordered models. 

3.1. Distance on probability measures. Let us first consider the set A4x of probability 
measures on C([0,T],^f) with finite moments of order k (where k ^ 2 is given in (|2.1ip ) 
and endow this set with the Wasserstein metric 



<f?(Pi,P2) :=M\E 



sup 

s < T 



(3.1) 



where the infimum in (|3.ip is considered over all couplings $( 2 )) with respective 

marginals p\ and p2- Here, the i?W are understood as random variables on a certain 
probability space (fi,P). Note however that the definition of (|3,ip does not depend on 
its particular choice. (|3.ip defines a complete metric on M.x encoding the topology of 
convergence in law with convergence of moments up to order k (see [351 Th. 6.9, p. 96]). 
We endow M. x with the corresponding Borel u-field. 

Fix some probability measure m on C([0, T], X) x £ x [- i, i] (endowed with its Borel 
cr-field) that is absolutely continuous w.r.t. fj,( du) (x) dx. Thanks to a usual disintegration 
result (see e.g. [12l Th. 10.2.2]) one can write m as 

m(d9, du, dx) = m ul ' x (d6)fi{duj) dx, 

where (u, x) >— ► m u,x { d6) is a measurable map from fx [ — ^, (endowed with its Borel 
(T-field) into Mx- We consider the set M. of such measures m such that for all (co,x), 



m 



belongs to Mx, endowed with the following metric: 
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Definition 3.1. Fix p to be equal to 2 in the case of Hypothesis \(H1 )\ or as in (|2.33p in 

the case of Hypothesis \(H2)\ Then defir 



me 



Vmi,m2 6.M, 5T(mi,m2) := 



(3.2) 



Note that, by construction (see (|2.21[) ). the initial condition dvo(9, uj, x) = £( d9)fi( duj) dx 
belongs to Ai. 

3.2. The nonlinear process. The proof of Proposition 12 . 81 is based on a Picard iteration 
in the space Ai endowed with the metric introduced in Definition 13. 11 For fixed uj e 8 and 
Brownian motion B in X , independent of the sequence (-£>&)& ^ i, and for a fixed m e Ai, 
consider the following stochastic differential equation in X: 

d8{t) = c(9(t),u;)dt + J T(9(t),u,e,u>)^(x,x)m t (de, du), dx) dt + a ■ dB(t), (3.3) 

with initial condition 8(0) ~ Note that the previous system is well-defined, thanks to 
the regularity properties of T and c and since 



I' 



T(9(t),uj,9 1 ui) I ^(x,x)m t (d9, du, dx) 



< lirL f ld f m^(d0>(d^) dx^ ir^^*), 

J -4.4 Wxf 



=i 

where the quantity 

5(*):=supf d ^(x,x)dx (3.4) 

is smaller than 1 in case of Hypothesis |(H1)| and smaller that I±(^) (using (|2.15p ) in 
the case of Hypothesis |(H2)| Note also that no particular assumption on the probability 
measure /i is needed here, since we work in (|3.3[) for a fixed uj. 

Let us denote by : Ai —> Ai the functional which maps any measure m(d8, duj, dx) e 
Ai to the law 0(m) of (9, uj, x) where (#t)o < t < t is the unique solution to ()3.3|) . Note that 
the functional effectively preserves the set Ai. Proposition 12.81 is a direct consequence 
of the following lemma: 

Lemma 3.2. The functional admits a fixed point v in Ai. 
Proof of Lemma VS. "A As in |33| , we prove the following 

Vrai,m 2 EM, Vt<T, 5 t (0(mi), 0(m 2 )) K < C T f S u (m 1 ,m 2 ) K du. (3.5) 

Jo 

If (13. 5p is proved, the proof of Proposition 12.81 will be finished since in that case, one can 
iterate this inequality and find 

Vk > 1, 5 T (e k+1 (u ),e k (u )r < C^8 T (Q(uo),u ) K , 

which gives that {® k ( l/ o)) k> , is a Cauchy sequence, and thus converges to some fixed- 
point v of B. Let us now prove (|3.5p . The key calculation is the following: there exists a 
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constant C > such that for all 9%, 02 e X, oj e £, x e [ — |, for all mi,m2 £ 



Indeed, 



J r w, ■) •) dmi,t - Jr (0 2 , w, -, •) *(x, •) dm 2 ,t 



<C(||0 2 -0i|| a l + 5((mi,m 2 )). (3.6) 



<Jr < J r (0i, a;, •) *(x, •) dmi,( - J r (0 2 , w, ■) * (s, •) dm M 
J r (02, u), •, •) •) droi, t - J r (0 2 , w, •, •) *(a?, •) dm 2ji 



(3.7) 



:= <yri + <*r 2 . 



The first term <5Pi in (13. 7f) is easily bounded by || T \\ Lip S(^) || 2 — 0i ||, where S(^) is 
defined by (13. 4j) . The second term ^2 in (I3.7jl can be successively bounded by 



<5r 2 



f *0M)( fr (02,0;, 0,u>) m£f(d0)- \v (9 2i u>, 

M x Jr(0 2 ,o;,0,a))m? j f(d0)-|r(02,a;, 



9,oo) m 2 ' t (d0) J dx/x(du;) 



0, u>) m^fi d0) 



dx/i( dw) 



Note that the first term in the last inequality is always bounded: it is straightforward in 
the P- nearest neighbor case and comes from Remark 12.161 in the power-law case. Indeed, 
q has been precisely chosen so that qa < d, so that ^(x, -) q is integrable. 

Using the Lipschitz-continuity of T, we see that, for any coupling m £J ' a: (di?i, d$2) of 
m 1 and m 2 

ST 2 < C || T \\ Lip M ^ (E m .,, II Mt) - Mt) \\) P dx^diDU " , 

< c II r L ip ( f t ([ew,. || tf x (t) - * 2 (t) fp) P dx/i( d*)V. 

By Definition 13. 1\ this gives <5r 2 ^ C \\ T \\ Lip St(rni, m 2 ), which proves (|3.6p . We are now 
in position to prove (|3.5|) . Let us consider (0i,w,x) and (9%,oj,x) solutions to (|3.3|) for 
two different measures m\ and m 2 in Ai driven by the same Brownian motion, with the 
same initial condition. We have for all ^ t ^ T, 



l(t)-02(t) 



Jo 



(s) , c(0i(s),w) - c(9 2 {s),u})) ds 



+ 2^ <^0i(s) -0 2 (s), Jr(0 1 (s),o;,-,-)^(x,-)dmi - Jr(0 



2(5), uj, -, •) ^(2, •) dm 2 ) ds. 
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Using the one-sided Lipschitz condition (|2,9|) and (|3.6p . we obtain 

\\e 1 {t)-9 2 {t) || 2 < C \ \\9 1 {s)-9 2 {s) || 2 a l&s + C \ \\e 1 (s)-e 2 (s) \\5 s (m x ,m 2 )ds, 
Jo Jo 

C [ \\9 1 (s)-9 2 (s) || 2 a Ids + C f 5 s (mi,m 2 ) 2 ds. 
Jo Jo 

Consequently, using GronwalFs Lemma 

sup \\9 1 (s)-9 2 (s) II 2 a 1 ^ Ce CT f 5 s (m 1 ,m 2 ) 2 ds. 

s Jo 

Elevating this inequality to the power § ^ 1 gives 

sup||0i(«)-0 2 (a)irAl< {Ce CT Y {{ 6 s (m u m 2 ) 2 dsY 

^ {Ce CT YT^ f 5,(roi,m 2 )"da, 
Jo 

which gives 

^ (6(^1)^,6(7712)"'*) ^ {Ce CT YT^ ([ d s (m 1 ,m 2 } K d s y . 

Elevating this inequality to the power p and integrating over u and x leads to the desired 
result (13. 5p . Lemma 13, 21 is proved. □ 



We are now in position to prove Proposition E 

Proof of Proposition \2.8[ It remains to prove that if v is a fixed point of O then v is a 
solution to the weak formulation of the continuous limit ()2.22p . Indeed if v = 6(z^), one 
can write v(d9, duj, dx) = i> L),x { d0)/i( duj) dx where, for fixed u,x, i? ul,x (d9) is the law 
of the process solution to (|3.3p . Applying Ito's formula, one obtains for all f(9,co,x), C 2 
w.r.t. 9 with bounded derivatives, 

/(0(f), u;, x) = /(0 O , u,x) + \ f div e (aa T V f) (9(s),lu, x) ds + f V e f ■ c(9(s),lu) ds 

1 Jo Jo 

+ f V e /- fr(0(t),w,0,w)^(x,x)^(d0)/x(da))dxds+ f V /(0(s), w, x) • (<r d£ s ). 
Jo J Jo 

(3.8) 

Integrating (|3.8p over z/(d0, dw, dx) leads to (I2.22p . But in order to do so, we need to 
know that the term Vef (0, w, x) • c(0, w) is integrable w.r.t. the measure j/ aJ > x '( d0)//( dw) dx 
(the other terms are integrable, by assumptions on /). This is ensured by (|2.1ip . the fact 
that (by construction) u U] ' x {d9) has finite moments up to order k, and the fact that \x has 
finite moment of order 1 (recall (I2.12P ). □ 

The rest of the document is devoted to provide a proof for Theorems 12.121 and 12.171 
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4. Definition and properties of the propagator 

For reasons that will be made clear in Remark 14.21 below, we make in this section, 
as well as in Sections and \6\ some supplementary assumption on the regularity on the 
dynamics c: 

Assumption 4.1 (Additional regularity on c). We assume that for all u, the function 
9 >— > c(6,u)) is globally Lispchitz continuous. 

Of course, the FitzHugh-Nagumo case does not enter into the framework of Assump- 
tion 2J3 The purpose of Section [7] will be to discard this assumption. 

In this section, the function is either defined as in (|2.13p or as in (|2. 14[) . We know 
from Proposition 12.81 that there exists at least one measure-valued solution t >— >• v t to 
the continuous equation (|2.22p . We fix once and for all one such solution. We can then 
consider the stochastic differential equation: 

d6(i) = c(6(t),oj)dt+ lr(e(t),u,e,ij)^(x,x)i/t(de, dw, dx)dt + a- dB(t) 
= : c(6(t),u)dt + v(t,9(t),uj,x)dt + a ■ dB(t), 

where 9(0) ~ £. Thanks to the regularity properties of T and c and to the integrability 
of "f, (|4.ip has a unique solution. Define the propagator corresponding to (14. lj) : 

Vs,te[0,T], P s>t f(9,cu,x):=E B f(& s (9;cu,x),oj,x), (4.2) 

where <3>* (9; oj, x) is the unique solution to (|4.ip and E# is the expectation w.r.t. the 
Brownian motion B. 

Remark 4.2. If f is C 2 w.r.t. the variable 9, under A s sumptions 1 2. 2\ and \4-l\ made about 
c and r, it is standard to see that the propagator P s tf is of class C 2 in 9 and C 1 in s and 
satisfies the Backward Kolmogorov equation (see for example [30j Remark 2.3]j: 

V(0,cj,x,s,t), d s P S}t f(9,LU,x) + ^div6> (aa T W e P Stt ) (9,u,x) 

+ {[c{9, u) + v{t, 9, u, x)] ■ V fl ) P s ,tf(0, oj, x) = 0. (4.3) 

The main problem which motivates the work of Section [7| at the end of this paper is that 
proving similar Kolmogorov when Assumption ^. 1\ is discarded appears to be difficult. Nev- 
ertheless, we work in this section under this additional hypothesis and we provide in Sec- 
tion [?] a way to bypass this technical difficulty. 



The key calculation of this work is the object of Lemma 14.31 
Lemma 4.3. Let f : X x £ x [ — |, i] d — * H, be C 2 w.r.t. the variable 9. Then, 

i r T 

f, 4 N) -v T ) = (P , T f, 4 N) ~ Vo / + \K^\T l ] Q V e (Pt,Tf)(9 k (t),oj k ,x k )io-dB k (t)) 
1 C T 

-rj-rY. Ve(Pt,Tf){e k (t),0J k ,X k )-\(r(9 k ,0J k ,;-)V(x k ,-) , v[ N) - V t \] dt. (4.4) 
l-M ^ Jo L\ /J 
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Proof of Lemma \4--^\ An application of Ito's formula gives: for all k and < t < T, 

Pt,rf (9 k (t),u) k ,x k ) = Po,rf (Ok(0),u} k ,x k ) + d s P S)T f (9 k (s),oj k ,x k ) ds 

Jo 

+ f V e P s>T /(6'fc(s),Wfc,x fc ) • d9 k (s) + \ [ div e (cra T \7 e P s ^f) (9 k (s),u k ,x k ) ds. 
Jo z JO 

Using the definition of 6 k (recall (12.2[) ) and (14, 3D we obtain: 

Pt,rf {9 k (t),u k ,x k ) = P Q)T f {9 k (0),oj k ,x k ) - v(s,9 k (s),u k ,x k ) ■ V e P s ,Tf {9 k {s),u k ,x k ) ds 

Jo 

+ J V e P s , T f (9 k {s),u; k ,x k ) ■ (T (9 k ,u> k , ■, •) f (x k , •) , ds 

VeP Stt f(9 k (s),uj k ,x k ) ■ (crdB k (s)). 
Then, using the definition of v (•) (recall (|4.ip ) and summing over k lead to: 
Pt,xf, 4 N) ) =(Po,rf, 4 N) )+ VePs,tf(9 k (s),oo k ,x k ) • (adB k (s)) 



f 

Jo 



-Ef 



it 



A straightforward calculation using (|4.3p shows that (PtTf , v t) = 0. Using this and the 
previous equality, one obtains the desired result (choose t = T and recall that Pt,t$ = /)• 
Lemma 14.31 is proved. □ 

The purpose of the following lemma is to establish regularity properties of the propa- 
gator P t ,rf- 

Lemma 4.4 (Estimates on the propagator Pt,T)- Fix T>0,0<t<T and a e [— i, 

(1) Assume VI/ satisfies Hypothesis \(H1 )\ For any R e (0, 1] and any f in C^ a , Pt,rf 
is also in C^ a and one has the following estimate 

II Pt,Tf\\ R>a < ^IH^II/II^, (4.5) 

/or some constant \\\P\\\ (that can be chosen equal to L + 3/2 || T \\ Lip , recall (12.9P ). 

(2) Assume ^ satisfies Hypothesis (H2)\ For every a e [ - i,^]* 1 , for any f in C a , 
Pt,rf is also in C a and one has the following estimate 

\\Pt,Tf\\ a <M\e m ™Ul, (4-6) 
for some constant \\\P\\\ (that only depends on T, ^ and c). 

Proof Note that, by a usual density argument, one only needs to prove (|4.5p and (|4.6|) 
for test functions / that are C 2 w.r.t. 9. Fix T>0, 0<i<T, ae[— and consider 
two different flows for (14. ip $ s (9i\ Ui, x), for i = 1,2, with different initial condition and 
parameter but at the same site x, with the same Brownian motion. For simplicity, we 
write instead of $> s (9i; ui{,x). Then, using the one-sided Lipschitz condition (12.91) on 
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c, we obtain 

I **(2) - || 2 = || 9 2 - 9 1 || 2 + 2 f <$«(2) - , c(<&?(2) )W2 ) - c(*«(l), Wl )> dn 

+ 2 f <*«(2) - , «(«, <^ (2), W2 , x) - v(u, Wl , z)> dn, 

s£ || 2 - 0i || 2 + 2L J (|| *?(2) - $«(1) || 2 + || w 2 - Wl || 2 ) dn 
+ 2 f || ^(2) - <D«(1) || || W ( U) $^(2),a;2,x) - <"l>*) II dn, 

:=Sv(u) 

where the definition of v(-) is given in (|4.ip . The Lipschitz-continuity of T implies 

< J I r($?(2), wz, 0, w) - r($*(i), wi , 0, a) || d0>( d ^> ds 

< II r H^.srC*) (|| *?(2) - || + |] a7 2 - w a ||) , 

where S(*&) has already been defined in (|3.4|) . Putting things together we see that, for 
C = 2L + 3||r|| Wp 5(*), 

I $*(2) - $*(1) f ^ ||02-0i|| 2 + c£(||^(2)-^(l)|| 2 + || W2 -a;i|| 2 ) dn. (4.7) 
An application of Gronwall's lemma leads to 

I $*(0 2 , W2, x) - & s (9i, U!,x) || 2 + || oj 2 - wi || 2 < e c(t - s) f || 9 2 - a || 2 + || oj 2 -u^f 



(4.8) 

Then, in the case where \P satisfies Hypothesis (Hl)| we have P^t/(0, u, x) = xr( x ~ 
a)g($J(Q;u),x),Qj), when f(0,u,x) = xr{x — a)g(9,uj). But then, 

I g{$J(9 2 ; U2,x),U2) ~ g(^f(9 1 ;u 11 ,x),u 11 ) f < || / ||^ a (|| $f(2) - ^(1) || + || a; 2 - ||) 2 , 

^2 || /|| 2 (l *n2)-*f(l) r + 1 || 2 



< 2 II / llfl,a eC(T ~ t} ( II ^ " 9 l || 2 + || ^ 2 " Wl || 2 

so that 

I g($T(e 2 \ U2,x),u 2 ) - g(^f{9 1 ;u 1 ,x),oj 1 ) || < V5 1| / || R>a (|| 2 - 0! || + || w 2 - wi 

which is the desired estimate (|2.24p and gives (|4.5p . The same kind of calculation in the 
case of Hypothesis |(H2)| leads to the estimate (|2.27p for Pt,rf- 

Thus, it remains to prove estimates (12.281) and (|2.29p for P* t/ i n the case of Hypoth- 
esis |(H2)| The case of (|2.28|) is straightforward. As far as (|2.29p is concerned, the same 
kind of calculation with two different flows & s {x) := & s (9;uj,x) and & s (y) : = $*(0;u;,n), 
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with the same 9 and uj but at different sites x and y leads to 
I $*(x) - || 2 < 2L f || *«(x) - C(y) II 2 dn 

J s 

+ 2 f || *«(x) - || || v(u, w, x) - r;(«, ^(y), w , y) \\ du, 

Js " v ' 



=8v(u,x,y) 



with, 



Sv(u, x, y) < J I r(^(x), w, 0", x) - r($?(y), w, 9, x) || d0)/i( dcD) dx 

< J I r(c^(x),a;, 0", w) - m u s (y),oj, 9, Q) || *(x, s)j# a (d0)/i( da;) dx 
+ J I m u s (y), uj, 9, uj) I x) - vl/ (y, x)|<- s ( d0)/i( da;) dx 

< lirii^s-Wdi^x)-^)!!) 



f |*(x,x)-*(y,x)| f C s (d0)Mda)) dx 
Jr-i,il d Ja"x^ 



=i 



id— a 



< || r || Lip sm || - <^(y) || + z 2 (*) || r || x - y 

where is defined in (|3,4p and where we used assumption (|2.16p . This gives, for 

C = 2L + 2||r|| Lip 5(vI/)+X 2 (M/)||r|| 00 , 



ft 

J i ,.i ;tJ i . ,i II 2 ^ n \ || ^tt/„, ;un , ... ii-' 



x) - $?(y) || 2 du + Z 2 (*) || T 11^ (t - s) || x - y 



|2(d-a) 



Consequently, by Gronwall's lemma, 

I $1(6; uj, x) - $1(9; uj, y)f ^ Z 2 (*) || T \\ x (t - s)e c ^ \\x-y || 2 ^ 
Then, for any < t ^ T, we have 

|| rf/ ',.-,•/ || 2 := J || x - a || 2 ^ D.rfiO. uj,x)-\\y-a ^ P tjT f(9, uj, y) f 

= | \\x-a\\ 2 ~< f(<f>J(9;uj,x),uj,x) - \\ y - a f" 1 f($J(9; u,y), uj, y) 

< f || x - a || 27 I f($J(9; uj, x),u, x) - f($J(9; uj, y),u>, x) \\ 



(4.9) 



x - a \\^ f($J(9; uj, y), w, x) - || y - a j| 2 ^ f($f(9; uj, y),u>, y) 



< \\f\\l(j$f(x)-$J( y ) 



+ || x - y 



127— a 



2(27-0) 



^2||/|| 2 (Z 2 W lini^CT-t) vl)e' 



C(T-t) 



x - y 



i2(d-a) 



(4.10) 



x - y 



2(27-a) 



(4.11) 
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where we used assumptions (|2.28|) and (|2.29|) in (|4.1U|) and the estimation (|4.9|) in (|4.11|) . 
Using the definition of 7 (recall (I2.18P ). it is always true that d— a ^ 27 — a. Consequently, 



x ~ a II 27 Pt,rf(0, u,x)-\\y-a || /7 Pt,rf(0, u, y) 



2"/ 



(T-t) 



f\\ a \\ x -y 



1 27— a 



which leads to (|2,29p . Lemma 14,41 is proved. 



□ 



Remark 4.5. One could wonder why we have not simply used in the calculation above the 
global Lipschitz assumption about c (recall Assumption \4-l\ ), instead of the more involved 
one-sided Lipschitz inequality used here. The crucial reason for this is that in order to be 
able to discard Assumption in Section [?] below, we need to ensure that the estimates 
of Lemma \4-4\ do not depend on the modulus of continuity of c, but only on its one-sided 
Lipschitz constant L. 

Using (|4.5p (respectively ()4.6|) ) in (|4.4p . we easily see that for every a e [ — i, for 
any given / e C^ a with \ f \ Ra ^ 1 (respectively / e C a with || / || a ^ 1), we have 



f,4 N) ) 



{f , v T ) II < II ( P ,t/ , ^ ) ) - (Po,rf , vo) 



+ 



J— Y{ T ^o^Tf)(O k (t),u; k ,x k )-(a dB k (t) ) 
l A ivl IT Jo 



1 r T 



V e P ijT / I 



(r(e k ,oj k ,;-)^(x k ,-) , 4 N) - u t y\\ at. (4.12) 



Using (|2T25|) and (USD (resp. ([2T3TJ) and ((Ml)), the term || V B Pt,Tf II (O k (t),u k , x k ) in the 
third summand of (|4.12p can be bounded by V2e'" p "^ T_ ^ ||xr||oq m case °f Hypothe- 
sis [(HI)] and by I Xk — a || a \\\P\\\e^ p ^ T ~ t ^ in case of Hypothesis |(H2)[ In both cases, the 
bound that we find can be written in the form 



V e Pt,r/ || (9k(t),u) k , x k ) ^ e 



\\P\\\(T-t) 



p(x k ) 



(4.13) 



(p is a constant in the first case and proportional to \\x k — a \\~ a in the second). In 
particular, it is uniform in / and (6 k ,LJ k ). Let us now fix the integer p equal to 2 in the 
case of Hypothesis |(H1)| or defined as in (I2.33P in the case of Hypothesis |(H2)[ Elevating 
inequality (|4.12p to the power p and taking the expectation lead to 



1 



3P-1 
+ E 



E 



E 



f,4 N) 



< E 



1 



|Avf? 



f Vg(Pt, T f)(0k(t),CJ k ,x k ) • (<rdB fc (t)) 
Jo 



1 Vj" ei p *Wp(z k )\(T(0 hl u, kl -,.)9( Xkl -) , v[ N) 



AjvI^. 



(4.14) 



At 
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Let us concentrate on the third term of the last inequality, that we denote by D^. By 
successive use of Holder's inequality (recall that | + | = 1), one has: 



D N ^ 



j; e «-) d4 )' E j;^ ? ^, K r ( , 



9 k ,LU k , ■, •) * (x k , •) , v\ 



(JV) 



dt 



3 V\\\P\\\T _ 1 

""MP 



E 

9 r T 



l A iv| V /Jo \ A N\ V 11 X 



Wfc, •) (xjfc, •) , V, 



(JV) 



(4.15) 



At this point, here are the main steps of proof that we will follow in the remaining of 
this paper: we have built the spaces of test functions (recall Definitions 12.91 and 12. 13[) in 
such a way that they precisely include the functions (9,uj,x) >— ► r (9k, Uk, 9, (J) \P (x&, x) 
for all h (in this case, a is equal to x^). Since the distances between two random measures 
introduced in Definitions 12.111 and 12.151 are exactly the suprema of evaluations over all 
such test functions, we are thus able to bound the term within the integral in (|4.15p in 
terms of the distance between and v. 

The second point of the proof is to obtain an estimate (uniform in /) of the speed 
of convergence to of the two first terms in (|4.14p . Taking the supremum over all test 
functions / and applying Gronwall's Lemma lead to the conclusion. 

Those steps are somehow easy to follow in the P-nearest neighbor case (see Section [5|) 
but are more technically demanding in the power-law case (see Section [6]). 



5. Law of Large Numbers in the P-nearest neighbor case. 



The purpose of this section is to prove Theorem 12.121 Thus, throughout this section, 
we suppose that ^ satisfies Hypothesis |(Hl)| for some R e (0, 1]. In this case, the integer 
p introduced in (|4.14p is equal to 2 and the function p in (|4.13p is bounded (equal to 
\f2 || xr IIqo)- Ln particular, the two terms in front of the integral in (|4.15p are trivially 
bounded by a constant, equal to 6 II XR 



2 

loo" 



The following proposition proves the convergence to of the first term in (I4.14p together 
with explicit rates: 



Proposition 5.1 (Convergence of the initial condition). There exists a numerical constant 
C\ > (independent of R) such that for all f e (J r -\dCn a with ||/||p a ^1 and 

II /L <i 



E 



(P , T f,V { N) )-<P0,Tf,V0) 



Ci 

J\[d/\2 ' 



(5.1) 



Proof of Proposition \5.1[ Recall that the couples (0j(O),u;j)i ^ j ^ at are supposed to be 
i.i.d. samples of the law £( d9) ®m( dw) on X x £. Let / e Cr^: by definition, f(9,uj,x) = 
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g(9, oj)xr{x — a) so that -Pq.t/ = x( x ~ a )Po,T9- Let write ip := Pq^tQ for simplicity. Then: 



S, (f ) , E <P ,t/, 4 N) )-(Po >T f, vq> 



E 



n-|E^(^ w iW x i - a) - I <p(6,ui)xr(x - a)C( d6»)/i( dw) dx 
2E XflCxj - (Wi^i) - J^,a;)C(dflHcL;) 

J ¥>(0, w)C( d0)/i( dw) ^l^i 2 XHfo ~ a ) ~ J ^(^ ~ a ) dx 

2 (vWjW) - J^^)C(de)Md^) 



(2i2) 



2 || 9? | 



An + P>n . 



(5.2) 



Since the (0j,cUj) are i.i.d. random variables (with law C® a standard calculation shows 



A 



N 



\A N \ 2 (2R)^ 



£e (p(e jtUj )- j (p(0,u)<Z(M)n(du) 



8||/ 



2 d 7V d ' 



since || <p L = || P ,T5 L = (2i?) d || / ^ and |Ajv| = (2JV + l) d > {2N) d . 

Let us now turn to the case of the term in ()5.2 j) . We place ourselves in the case of 
non-periodic boundary condition (recall § I2.2.4p . The periodic case is simpler and left to 
the reader. Let a = (ai, . . . , cy)- One has 



>[-M] 

In the same way. 



1 \x-a l \ sSfldx 



(5.3) 



1=1 



1 d / 1 N \ d 

ia^s^(^ - «) = n ^ 2iz( 2 jv+i) j jy i '^- a '' < *J := n^^)- 

Then, from the obvious equality 

d d d 

[JZjv(aj) - jJZ(aj) = ^] Zjv(ai) . . .Z N {a k -i) (Zjv(a fc ) - Z(a fe )) Z(a fc+ i) . ..l{a d ). 



i=i 



i=i 



k=l 



and a recursion argument, one only needs to consider the case d = 1 in order to prove 
(15. ip . An easy calculation shows the following: for all a e [— ^ , 2], for all R e (0, l], 

^ (P + ± + a) , if-£<a<-£ + -R, 



1 fi 2iJ 
Z(a) = — J ^ h dx = j 1 



2R 



(R+l 



if -i + P^a^i-i?, 
if |-P< a< |. 
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Thus, in the one-dimensional case, we need to distinguish three cases, depending on the 
position of a e [— \ , 5] w.r.t. i?; we only treat the case — ^ — \ + R, the two others 

being similar and left to the reader. In this case, one has successively, 



\l N {a) -1(a) 



1 



AR 2 
1 

4R2 



1 N 
2N + 1 .^-1 



j-2oJV| < 2fliV 



1 



2iV + 1 



(i? + a) 2 



([2iY(fl + a)J +JV) - 
< (2^-1/2)2 < 1 



i? 



1 

" 2 
1 

2 " 



4i? 2 (l + 27V) 2 16R 2 N 2 
Proposition 15.11 is proved. 

We are now in position to prove Theorem 12.121 



AN 2 ' 



□ 



Proof of Theorem \2.1SX Fix some a e [— ~ , I] and some / e Cr jQ such that || / \\ Ra ^ 1 
and || / || oo ^ 1. Let us first give an estimate of the second term in (|4,12j) . Recall that B k 
is a Brownian motion in X = R m so that B k may be written as m i.i.d. Brownian motions 
(i) B H 

k ' • • • ' D k 
Lemma 14.41) in (15.5 



5,^V • • , 5,Y" ; ) . Then, using ([235]) (recall Remark [HUD in ([53]) and using (03]) (recall 



E 



1 f T 



V e (P t)T /) {6 k {t),u k ,x k ) ■ &B k {t) 



1 ™ ~7 ■ 
|Ajv| fc ,ti Jo 

^^^ril^T/llLdt (5.4) 
Jo 



|A 



N 



>n || XR 



IA 



i2j T e 2 lll p IIK T -')dt (5.5) 

m( e W-i) c 2 
(212) M |A^| ^ iV d ' { } 

where C2 = m<ye %d R \d~ 1 ^ an d where |||P||| is defined by (|4.5j) . 

Let us now give an estimate of the term Djy in (|4.15j) : by Definition 12.91 due to the as- 
sumptions made on T, it is easy to see that for fixed k the function f k := F (9 k , ou k , •) \I/ (x^, •) 
belongs to Cr,^ with norm || || J j a . fe = || F \\ Lip . Consequently, by construction of the dis- 
tance dn (recall Definition I2.11h . one has the following: 



Vt > 0, E 



(F(e k ,oj k ,-,-)^(x k ,-) , u[ N) -u t ) 
Putting together (|4.14j) . (|5.ip and (|5.6p . we obtain finally 

n C 2 



E 



f,4 N) 



2 Ci 
«S 3 



eW-1 2 
3 (2P) M |||P||| l|r " Li P 



Jo 



^) 2 dt. 



Taking the supremum over all functions / in Uoe[-i i] d ^R> a anc ^ a PPlyi n g Gronwall's 
lemma leads to the result. Theorem 12. 121 is proved. □ 
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6. Law of Large Numbers in the power-law case 



We suppose in this section that the weight ^ satisfies Hypothesis |(H2)[ 

Let us begin with a technical lemma that will be of constant use throughout this part: 

Lemma 6.1. There exists a constant Cq > (that only depends on j3), such that for all 
N,K ^ 1, for all a e V K , 

(1) for all < /3 < d, one has 



2 



J 

2N 



<c {N d K d ifa$V N 
^ °[N d ifaeV N . 



(6.1) 



(2) for P = d, one has 



2 

j; j/N^a 



J 

2N 



-d 



^C 



K d N d In N ifa$V 
N d hxN 



N 



if ae V N . 



(6.2) 



(3) for all (3 > d, one has 



y — 

V 2N 

j; j/N^a 



~ P <(J {N^K? ifatV N 
^ °\nP ifaeV N . 



(6.3) 



Remark 6.2. The estimates given in Lemma lfi.il in the case a e T>n are standard and 
optimal. The main technical problem of Lemma [6A\ lies in the case of a $ T>jy: in this 
case, the point a of the discretization T>k can be arbitrarily close to one point in the 
above sum. Those points belong to the discretization T>n. The minimal distance between 
a and the discretization T>n depends on K (actually it depends on the greatest common 
divisor of K and N, see the proof of Lemma \6.1\) . This explains the dependence in K of 
the estimations of Lemma \6.1\ 

The proof of Lemma \6.1\ is postponed to the appendix. Lemma \6.1\ will be at the basis 
of most of the estimations in this section. 

Theorem 12.171 is a consequence of the two following propositions: 



Proposition 6.3. Let fix a, 7 and p as in § \2.3l p. 0' there exists a constant C\ > 
(that only depends on p and Cq defined in Lemma \6. 1\) such that for all K ^ 1, N ^ 1, 
a e T>k and f e C a with \\ f \\ a ^ 1, 



E 



(P ,Tf,4 N) )-(P0,Tf,^) 



\Nf J 
K d lnN\ P 



d 



3d \P 

K~T In TV 



ifae [0,f), 



d 
2' 



1W 1 



if a 
ifae . 



(6.4) 



Moreover, in the case where a e T>n, the previous estimates are true for K = 1. 
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Proposition 6.4. Let a, 7 and p as in § \2.'A p. 0. There exists a constant C2 > sac/i 
i/iai /or all K ^ 1, for all a e for all f e C a such that \\ f \\ a ^ 1 



E 



|A 



( K±\z 
\N d ) ' 



*/ae[0,f), 
^C 2 U^)\ ifa = l 

(6-5) 

Moreover, in the particular case where a e T>n, the previous estimates are true for K = 1. 

Let us admit for a moment Propositions 16.31 and 16.41 Then the result of Theorem 12.171 
is a straightforward consequence of the following proposition: 

Proposition 6.5. Under the assumptions made above, there exist constants C3 and C4 
such that for all K,N ^ 1, one has: 



sup d K {v[ N \ vt) < C3 < 



i/ae [0,|), 



^ t sg T 



lnjV K d c CiK 



if a = 4, 



(6.6) 



where q in (j6.6|) zs iae conjugate of p and where C3 and C4 are large enough constants that 
depend only on p, T, V, ^> , c and on the constants C\ and C2 defined in Propositions RT 
and\6.4\ 



Proof of Proposition [675[ Let us fix K ^ 1, a e and / e C a with || / || a ^ 1. Let us 
recall the estimate obtained in (|4.14p and (|4. 15p : 



E 

3 p - x E 



(TV) \ P 



Po,Tf,v ( N) ~ v ) P +3P" 1 E ^2 V e (P t , T /) (9 k (t),u k ,x k ) • ((rdB fc (t)) 
7 |Ajv| -V Jo 



3 P 



-1 



3 9iii^r _ 1 



IAat 



. . r Zj ~, _ 



x k - a 



qot 



r T 1 

J j^2 E |( r (^' Wfc '-'-)^(^'-) >^ (JV) -^ 



dt. (6.7) 



We understand here the necessity of choosing p (and its conjugate q) different from 2. 
Indeed, the integer q (recall Remark 12. 16j) has been precisely chosen such that qa < d 

which ensures that the term Efc is «„i te; more precisely, an application 

dp 

of Lemma 16. 1\ (|6.ip shows that this quantity is smaller than K 1 whenever a e T>k and 
smaller than 1 in the particular case where a e T>n- 

Let us now prove (16. 6p in the case where K > N. Notice first that, thanks to 
the assumptions made on and T in § 12.31 for all k the function f k : (8,co,x) >—* 
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T (9 k , oj k , 9, <J) \P (x k , x) belongs to the space C Xk where x k e V^. Indeed (recall the defi- 
nition of Xi(*) (|2.15|) ). for all k and (9,u,9,6j,x), 

|| x - x k || a * (a*, x) I r(0 fc , w fc , 0, w) - r(0 fc , w fc , 0, w) I < Xi(*) || r || Lip (|| - I + 1| o - u I) , 
and 

|| x - x k \\ a ^ (x k ,x) \\T(9 k ,uj k ,9,uj) || ^Xi(^) I T 11^. 
As far as condition (|2.29p is concerned, we have (using (|2.17p ): 



x-x k f 7 f h (9,u,x) -\\y-x k f k (9,uj,y) 



2T 



^ II r 



\x - x k || 27 ®(x k ,x) -\\y- x k f 7 ^(x k ,y) ^ X 3 (*) || r 



\x-y 



27— Cf 



(p)/ 



Therefore, since K > N, by definition of the distance d K 
all fc, the following holds 



(recall Definition I2.15D . for 



E 



(T(O k ,U k ,;-)*(x k ,-) , v[ N) -V t ^J ^ m d%\u { t N \u t 



for the constant 771 := max ^Xi(^) || T || Lip , X\{^), T%{fy) || r . Using this estimate 

(|6.7p and taking the supremum over all functions / in (J ae x) L C a , one obtains 

i< L < K 



m 



4 ) (4' V) ^t) p < S^ 1 sup E I <P ,t/ , - (P 0>T f , vo) 

J—VCve (Pt, T f) (9 k (t),u k ,x k ) ■ (a dB k (t)) 

\ A N\ ^ Jo 

Jo 



3 p_1 supE 

/ 

-1 dp 



for r/2 := r\\ (^~^|pf^") * • The results of Propositions 16.31 and 16.41 together with an appli- 
cation of GronwalFs lemma leads to the estimate (]6.6p in the case where K > N. Note 

p-i 1 

that one can choose in this case the constants C3 := 3 * (2 max (C\, C2)) p (where C\ and 



C2 come from Propositions 16.31 and I6.4p and C4 : = 



3P- 



-Tr]2- 



Let us now turn to the case where K ^ N. In this situation, we cannot use Gronwall's 
inequality in order to obtain an analogous estimate on djr {i/( N \ v) , since the function 
f k (k e Ajy) defined at the beginning of this proof has not the sufficient regularity (f k 
belongs to C Xk where x k e T>n and hence may not belong to (J a eV K , C a for K < 



K 

1<K'<K 



N). Nonetheless, one can bound the term ^-E ( T (9 k , cv k , •, •) (x k , ■) , v£ — v t 



X N ) 



by supjE , vj. N ^\ — (/ , i/ t y , where the supremum is taken over functions / in 
UaeX> Ca with || / || a ^ 1. Using this estimate in (|6.7p and a calculation similar to the 
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previous one gives the following estimate: 



'(f)', if « 6 [0,|), 



sup sup 

0<t<T fe [j aeT , N Ca 



(e i </,^> -</,^>n < (c 3 c^) p u^f) 



d 

i ^p 



if a = |, 



(6.8) 

But then, for instance in the case a e [0, |) (we let the two other cases to the reader), for 
all K ^ N, for all / e (J ae x> K , C a for K' ^ K, inserting directly (|6.8j) into (|6.7p and using 
again Propositions 16.31 and 16.41 leads to: 



E 



</,,r>>-<>.->r^(£) P 



3' J 



-i 



q\\\P\\\ 



d \ V 



N'T 



Up to a change in the constant C 3 , this term is anyway smaller than (^K^J . 

Taking the supremum over all / in {J a eV K , C a , one obtains the result. □ 

K' < K 



The rest of this part is devoted to prove Propositions 16.31 and [ST 



Proof of Proposition \6.3[ Recall that the couples (#j(0), ^ j ^ jy are supposed to be 

Write again 



chosen i.i.d. according to the law C(d#) ® /i(dw) on x £. Fix a = -4 e f e C a 



with || / || a ^ 1 as well as a e (0, d) and the integer p ^ 2 as defined in 
99 := Pot/ for simplicity. Then, 

:= E J (p ,t/ , ^ } ) - <Po,t/ , ^o> |[ 



< 2 p ~ i E 



2 p 



-1 



v 

1 1 f 

2 J w ' x i^{ d9)n( duj) - J <p(6, uj,x)C( d0)n( du) dx 



'■= An + Bn- 

For simplicity, let us write Xj := ip(9j,Uj,Xj)— $ (p(9,u, Xj)((d9)/j,(dto); note that EXj = 
for all j. Since the (9i,coi) are i.i.d. random variables with law C ® A*> the first term An 
becomes 



A 



N 



IA 



, b/2j 



v 



2 2[p/2j 



i=l (fci + ... + fc; = [p/2j) ji,...Ji 
LP/2J 



2fy 

ii 



lA-JP Z_l ii _ i|2afci 

I N \ 1=1 (fc 1 + ...+fc i = [p/2j)ji,...,j i II x ii a H 



I 2ctfc; 



(6.9) 
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where we used || / || a ^ 1 and assumption (|2.28|) in (|6.9|) . Let us concentrate on the 
contribution of / = 1 to the sum in (j6.9|) . that we call An (where p = 2[p/2\): 



A N 



2 p 



| A ..ip 



AT ~ || Xj — a 



2pa ' 



Here, one has to distinguish two cases, depending on the value of a e [0, d): 

(1) If ^ a < i then by definition p = 2 and pa < d so that an application of 
Lemma [67T1 (|6.ip leads to 



(6.10) 



(2) If a ^ §, then p is chosen such that p > so that pa > d. Then Lemma 16,1 
leads to: 



i ivrp Q 

JV ^ ]\[pd u u ]yp(d-a) 



(6.11) 



It is also easy to see that the other terms in (|6.9p are negligible w.r.t. An as 
N -> 00. 

1 

Let us now turn to the second term Bn- {Bn) p is the difference between the Riemann 
sum of the function <3? := x 1— * \ (p(9,u;, x)£( d0)/i( dw) and its integral, so that it should 
be small with N. But one has to be careful since (p as a discontinuity (</? belongs to some 
C a for some a) and since we want to have a result uniform in the function cp: 



1 



2p-i 



iV 



p 

|^— I ^ $(xj) - J<I>(x) dx ^ XlJ || ^(xj) - $(x) || dx 



(6.12) 



where Aj := iz e [ — ^ , VA; = 1, . . . , d, jk ^ Zk < jk + ^ j ^ s the infinitesimal sub- 
domain of A^r of size 277 of corner j. Let us begin with the following straightforward 
inequality: 

|| $(x) - $(y) || < (I || x - a |p 7 



\y-a\ 



x - a || 7 <£(x) + || y - a || 7 $(y) 



1 



17 II 117 

x — a || ' || y — a \\ ' 



<3?(x) || x — a || 27 — $(y) || y — a 



l 2 1 



(6.13) 



Using the assumptions made on /, we deduce in particular from (12.28P and || / || a ^1 that 
|| x — a || 7 <I>(x) is bounded by || x — a || 7 ~ a . Using also (|2.29p . it is then immediate to see 
that 



$(x) - $(y) || ^ 



x - y\ 



117 II 117 

x — a || ' || y — a \\ ' 



x — a 



\x - y 



17— a 



\x-y\ 



y -a 



\n-a\ 



\x-y 



\2j—a 



17 II 17 ' 

x — a || ' || y — a \\ ' 



\x-y 



127— a 



II O! II 117 

\x — a || || y — a || ' 



17 II II O! 

x — a || ' || y — a \\ 



I 117 II 117 ' 

x — a || ' || y — a \\ ' 



(6.14) 
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Using (|6.14|) in (|6.12p . one obtains that 

B N s$ 2 p 




(6.15) 



3 

The first of the three sums in (|6.15p can be bounded by the following quantity: 



^ X • n\ ii" ii \\ct\ ii n7 II x x j P 

t-r 1 mm (|| Xj-\ — a \\ , || Xj — a \\ ) \\ Xj — a \\ ' 

—y 



j V II J A 1 7 II J II I II J II "^j 

1 

_/y d +7 ^-i min (|| Xj-i — a \\ a , \\ Xj — a \\ a ) \\ Xj — a || 7 

Let us once again distinguish three cases, depending on the value of a: 

(1) if ae [0, |), then a + 7 < d (recall (j2.18fl ). so that an application of Lemma |6.1| 
dSIU) leads to 

'AT7 



s n < C —. (6.16) 



(2) if a = \ , then a + 7 = d (recall (I2TT8I) ). so that Lemma EH (JOJ) gives 



2 

K d lnA 



4 1} < Co^^. (6.17) 



(3) if a e (f,d) then a + 7 > d, so that Lemma 16. 1( (|6.3p gives 

' jV<i-<* ^ °iV^ 



S n ^ Co-^- < Co-,-. (6.18) 



/OA 

The same calculation leads to the same estimates for the second term in (16. 15j) . A 

(3) 

very similar calculation also leads to the following estimate for the last term S N : 



S§ ] ^C {" " ~ (6.19) 

^ ifae[f,d). 

Combining estimations (|6.19p and (|6.10p (resp. (|6.1ip ) and (|6.16p . (resp. (|6.17p or (|6.18p ) 
leads to the desired estimation (|6.4p . The proof of the case where a e T>n is analogous 
and uses the estimates for a e T>n in Lemma IBTTl Proposition 16.31 is proved. □ 

It remains to prove Proposition 16.41 whose purpose is to control the martingale term in 
gj2D: 



jv^-a if a e [0, 2;, 



Proof of Proposition \6.4\ Fix some K ^ 1, a e T>k and /eC a such that || / || a ^ 1. The 
martingale M t N : = jx^j Z fc ^ V <? (P^t/) {Qk(t),u k ,x k ) ■ dB k (t) may be written as = 

jX^ZkZZl ftfym ( P t,Tf) (e k (t),u k ,x k ) dB®(t), where for all k, B k = (b£\. • • ,£?£ m) ). 
Consequently, its quadratic variation process is given by 

1 m r-T 

(M N ) t = — — V V \\d em p t , T f(e k (t),u k ,x k )fdt. 

\A N \ k i=1 Jo 
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Applying Remark 12.141 and Lemma 14.41 we have almost surely that 

<n^^f e2||p,(r - ()d * 

|Ajv| £ \Xk ~ a\ Jo 

An argument repeatedly used in this work shows that one can bound the quadratic vari- 
ation by Cjjj (respectively C K ^ N and C N ^ {d 1 a ) ) when a. < | (respectively a = | 
and a > |), for some constant C > 0. Then, Burkholder-Davis-Gundy inequality 
E (I M t N \\ p ) ^ C P E ((M^j) gives the result. Proposition El is proved. □ 

7. The case of a locally Lipschitz dynamics c(-) 

One of the key arguments of the proofs of Theorems 12.121 and 12.171 is the fact that 
one can derive a Kolmogorov equation (recall (|4.3p ) for the propagator P s tf defined in 
(14. 2p . Under Assumption 12.21 on the dynamics c(-) (one-sided Lipschitz condition and 
absence of global Lispchitz continuity), deriving such a Kolmogorov equation appears to 
be problematic (see in particular |20} I18|) . Even if such a result existed, we could not find 
a proper reference in the literature. 

One can bypass this technical difficulty and prove nevertheless Theorem 12.121 and 12.171 
by an approximation argument. We will suppose throughout this section that c satisfies 
only Assumption 12.21 

7.1. Yosida approximation. Let us denote for all (9,uj), c{9,ui) := c(9,u) — L9, where 
we recall that L is the constant appearing in the one-sided Lipschitz continuity assumption 
(|2H . In terms of c, §2JB reads: 

V(0,w), (S,u), (9-9, c{9,lo) - c{9,Cu)) ^ 0, (7.1) 

and, for example, the mean field evolution (|4.ip reads: 

d9(t) = c(9(i),oj) dt + v(t,9(t),u),x) dt + a ■ dB(t), (7.2) 

where v(t,9(t),oj,x) := v(t,9(t),u,x) + L9{t). 

For all A > 0, consider c\ the Yosida approximation of c (see [U App. A] for a review 
of the basic properties of Yosida approximations): 

V(#,w), c x {B,u) := c(Rx(\9),u), (7.3) 

for 

V(0,w), Rx(9,u) := (A - cCw))- 1 (9). (7.4) 

Consider now the solution 9\ of the following SDE (with the same initial condition and 
driven by the same Brownian motion B as in (17. 2|) ) 

d9 x {t) = 5 x {9 x {t),oj) dt + v{t,9 x {t),oo,x) dt + a ■ dB(t), (7.5) 

that is, the analog of (|7.2p where c has been replaced by its Yosida approximation. Note 
that one can proceed exactly in the same way for the microscopic system (|2.2p . From now 
on, whatever X may be, the subscript notation X x will refer to the analog of X when the 
dynamics has been replaced by its Yosida approximation. Note that we will most of the 
time drop the dependencies of the functions in u, for simplicity of notations. 

It is easy to see that if c and c x have the same regularity in 9 (see e.g. [8]). Moreover, 
c x has the supplementary property to be uniformly Lipschitz continuous. In other words, 
c x satisfies Assumption 12.21 as well as Assumption 14. H so that everything that has been 
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done before is applicable: Theorems 12.121 and 12.171 are true in the case of an interaction 
ruled by cy. 

sup dUf,v t , x ) ^CN-P, (7.6) 
te[0,T] v ' ' 

for d either equal to dji(-,-) or d£> (•, •) and /3 one of the appropriate exponent appearing 
in the formulation of Theorems 12.121 and 12.171 Note that the constant C in (I7.6D does 
not depend on A. Indeed, the assumption made in Section H] about the global Lipschitz 
continuity of c was made only to ensure the existence of the Kolmogorov equation. In 
particular, the modulus of continuity of c did not enter into the calculation made in 
Section [U the only dependence in the dynamics c was in its one sided-Lipschitz constant 
L (recall Lemma l4.4p . which is conserved by the Yosida approximation. In other words, 
every constant estimates made upon evolution (17.5P is independent on A. 

Now, Theorems 12.121 and 12 . 1 71 in our general framework are an easy consequence of the 
triangular inequality and the proposition: 



d(v£>,vM)^0 (7.7) 



Proposition 7.1. For all N ^ 1, 

SUP U I K, A . K, j 

te[0,T] v 7 

sup d(i/ t> \,v t ) ->A— oo 0. (7.8) 

te[0,T] 

The rest of this section is to devoted to prove Proposition 17.11 Let us begin with some 
a priori estimate: 

Lemma 7.2. We have the following a priori estimates 

supE( sup || x (t) || 2 ) < oo. (7.9) 
A>0 Ue[0,T] J 

and, 



PjsupJ || c x {6x{s)) f ds < oo j - I. (7.10) 



Proof of Lemma \7.S\ Let us first prove the first estimate (17. 9p : applying Ito formula, 

\\e x (t) || 2 = ||0a(O) || 2 + 2 f (8x(s) , cx(ex(s)) + v(s,9 x (s),oj,x)) ds 

Jo 

+ 2 f (6 x (s), dB(s)) + tr(aa T )t, 
Jo 

^ || 9 X (0) || 2 + 2 (|| 5(0) || + L + || r 11^ 5(*)) f || 9 x (s) || 2 ds 

Jo 

+ 2 f (9 x (s), dB(s)) + tr(aa T )T. 
Jo 
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Taking expectations and using Burkholder-Davis-Gundy inequality, we obtain that for 
some constant C > (independent of A), 

E ^sup || 9 x (s) f \ ^ E (j| 0(0) || 2 ) + tr(aa T )T + 2cJ E (^sup || A (u) f^j ds 

i \ 



+ 6tr(aa T ) 1 / 2 B^ || 9 x (u) f du) 

E (|| 0(0) || 2 ) + tr(aa T )T + 2C f E ( sup || 0. 

^ ' JO \u^s 

+ 18tr(aa T )T + -E I sup || 9\(u) || 
2 V < t 



x(«) II 2 ) ds 



which implies 

E f sup || A (s) H 2 ] ^2 (E (j 0(0) || 2 ) + mr(aa T )T^ + AC f E f sup || A (u) I) 2 ") ds 

and Gronwall lemma leads to the result. 

Let us now turn to the second estimate (|7.10p : define Y x (t) := 6 x [t) — a ■ B{t). Then, 
Y\ satisfies: 

dY x (t) = (c x (Y x (t) + B(t),u) + v(t,Y x (t) + B(t),u,x)) dt. (7.11) 

Clearly, 

|| Y x (t) || 2 = || Y x (0) || 2 + 2 (\y x (s) , c x (Y x (s) + a ■ B(s))) ds 

Jo 

+ 2 f (Y x (s), v(s,Y x (s) +a-B(s)), u,x) ds 
Jo 

< || y A (o) || 2 + 2 (|| 5(0) || + l + 1 r iu s(*)) f || y A ( s ) || 2 ds 

Jo 

+ 2 f <T A (s) , c A (cr • B(s))) ds 

JO 

^ || y A (0) || 2 + 2 (j| c(0) I + L + || r || S(tf) + J || c x (a ■ B(s)) f ds^j J || F A (s) || 2 ds 
taking the supremum in A and using 1^(0) = A (O) = 0(0), we have 

sup || Y x (t) || 2 ^ || 0(0) || 2 + 2 I C + f \\c x (a- B(s)) f ds] f sup || Y x (s) f ds 

A V Jo /JO A 

sC || 0(0) || 2 + 2 (c + J c(a • B(s)) || 2 ds^j J sup || Y x (s) f ds 
where we used the pointwise estimate || c A (0) || ^ || c(0) ||. Gronwall lemma gives 
sup || Y x (t) || 2 ^ || 0(0) || 2 exp (2 (c + J || c{a ■ B(s)) f ds^j T^j 
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that is almost surely finite, since c is locally bounded and the trajectories of B are almost 
surely bounded. Consequently 

sup sup || 6\{t) || 2 ^ sup sup |[ Y\(t) || 2 + sup || B{t) \\ 2 < oo, a.s. 
Since c is polynomially bounded, this implies now that 

sup || c\(6\(t)) || 2 dt < oo, a.s. 

A Jo 

which is the result. □ 

The key estimate of this section is the following 
Proposition 7.3. Almost surely, the following holds 

limsup sup || 9(t) - 9 x (t) || = 0. (7.12) 

A^oo te[0,T] 

Proof of Proposition \7.3\ Let us fix A < [i. Since the Brownian motion is the same, one 
has successively (for a constant C = L + \\T \\ Li S(^)) 



■±e- 2Ct || 9,{t) - 9 x {t) || 2 = -2Ce- 2Ct || 0„(t) - 9 x (t) f 



+ 2e- 2Ct (9^t) - 9 x {t) , c^O^t)) - c x (6 x (t))) 
+ 2e- 2Ct (9^{t) -6 x (t), v(t,9^t),Lu,x) - v(t,9 x {t),u,x)) 

< -2Ce- 2Ct \\9 ll (t)-9 x (t)\\ 2 
+ 2e- 2Ct (9 fl (t) - 9 x {t) , c M (fy(t)) - c A (0 A (i))> 

+ 2e- 2Ct (L+\\r\\ Lip s(*)) II e,(t)-9 x {t) || 2 

< 2e- 2Ct <0 M (t) - 9 x {t) , c^t)) - c x {6 x {t))) 
= 2e- 2Ct (^ (RpW^t)) - ~c(i^(^(i))) 

R x (X9 x (t)) - ^c(R x (X9 x (t))) S j ,c(^(^(t))) - c(R x (X9 x (t))) 

< - 2e- 2Ct (h,{9^t)) - \c x {9 x {t)) , c„(0„(i)) - c x {9 x {t)) 



Integrating this inequality gives (since the initial condition is the same) 

(fy - 9 X ){T) || 2 < e-^Q-c^^t)) - ~c x (9 x (t)) , c M (^(t)) - c x (9 x (t))^ dt. 



1 —2CT " ' " " - 1,9 - ^ —ir<t / 1 _ / „ , .n% 1 

This gives in particular that 



J o e- 2Ct (jfMt)) - ^c A (0 A (t)) , ^(W) " cx(O x {t))^ dt < 0. 
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Let us denote as || • ||^ the Hilbert norm in H := L 2 ([0,T],e 2Cs ds; X). Then, from the 
identity 

i(cM - c x (d x ) , -c^) - jc x (e x )^ = Q + ~) || c^Op) - c x (e x ) f H 

one obtains that 

1 1\ - 2 / (I 1\ (\\ ~ fa m,2 „ s , fl m,2 



- + ^J ||5 M (^)-c A (0 A )||^ < ( X - - J (II £A) 11^-11^(^)11^) (7.13) 



which gives in particular that A h-> || c a (0 a ) ||# is increasing and by (|7.10p bounded and thus, 

~ C A (6>a) ||# ->A,A*-»co 0, 



convergent. The same inequality (|7.13p shows also that || c^{9^) — c x {9 x ) \\ 2 H — >a,^oc 0, so 



that (c x (9 x )(t)) converges in H to some c^{t). 

Going back to the first inequality of the proof, one has 

-2Ct /, 



I sup e- 2Ct |^(t)-^ A (t) || 2 < f e- 20t <e Ai (t)-W A (/).r > (^(/))-r A (# A (/))>d/ 
^ iero,Tl Jo 



+ T [ T e- 2Ct ||c M (^(t))-CA(0A(i))|| 2 dt 
Jo 

< ~ sup e- 2Ct || 0„(t) - A (t) || 2 
4 ie[0,T] 

+ T T e- 2Ct || c^(^(t)) - c A (0 A (i)) || 2 dt. 
Jo 



Hence 

sup e - 2Ci || ^(t) - e A (t) || 2 < 4T f T e - 2C * || c^e^t)) - c x (e x (t)) || 2 dt, 

te[0,T] Jo 

which goes to as A, fi — > oo. This implies that there exists an adapted process 9 with con- 
tinuous trajectories such that lim^oo 9 X = Q> uniformly and almost surely. Clearly, for all 
t, the strong continuity liniA^oc R\(^G(t)) = 9{t) of the resolvent and the uniform Lipschitz 
continuity || R x {\9 x {t)) - R x (X9(t)) || < || 9 x (t) - 9{t) || implies that lim A R\(\9 x (t)) = 
0(f). Finally, continuity of c gives lim A ->oo c x (9 x (t)) = c(R x (X9 x (t))) = c(9{t)). Conse- 
quently, we have that, almost surely c{9 t ) = (ko(t), so that 9 solves equation (|7.2p . so that 
by uniqueness 9 = 9 almost surely. □ 



We are now in position to prove Proposition 17.11 

Proof of Proposition 7. 1 . We only prove (17. 8h , the proof of (17. 7h follows from analogous 
estimates with the microscopic equation (|2.2p . We only treat the (more complicated) case 
of the power law interaction. Fix any / in C a for some a with || / || a ^ 1. Then, by 
Lispchitz continuity of / in the variable 9 

Kf , " </ . ^>l < W E £ || 9 x (t) - 9(t) || . 
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Taking the supremum in / and in t leads to 



sup d(v t , x ,vt)<S(y)E B sup \\6 x (t)-9(t) 

te[0,T] te[0,T] 



By (|7.12|) we have the almost sure convergence to of sup t6 r T i j| 9\(t) — 9(t) || and (|7.9[) 
gives the boundedness in I? implying uniform integr ability. The result follows. □ 



Appendix A. Proof of Lemma 16.11 
Proof of Lemma \6.1[ Let us proceed by induction on the dimension d. Let fix d = 1: 



Let us begin with the case where a $ V^: let J be the integer such that 



4rjtr- Then, an easy comparison with integrals shows the following 



2N 



< a < 



J 

2N 



-P 



J 



s^A^J \2aN 

rN 

+ \t-2aN\-P dt 



t\-P dt 



\2aN - J 



-P 



\2aN- (J+ 1) 



J+i 

2 /3 N /3 f \2aN 
Jo 

J 

2N 



t 



-P 



dt 



r N 

2 P N P 

J. J+1 



t-aN^ 13 dt 



J+1 



2N 



It is straightforward to see that the two first integral terms are smaller than 
whereas each of the two remaining terms is smaller than p(N, K)~@ , where 



p(N,K) : = 
and /3 < 1, 



inf 



\l\ < K 
j/N^l/K 



J_ 

2N 



I 

2K 



gcd(K,N) > 1 
2KN 2KN 



. Consequently, since K ^ 1 



E 



J_ 

N 



2N 
d-p 



2K /3 N f} < C NK. 



• The case where a e T>n is easier: in that case, a 
again by comparison with integrals, 

-P 



2^f for some k. Then, once 



J 
2N 



2V 2 \j 



k\ 



-P 



1-/3 



((/V + k) l - p + (N - k) 1 -^ 



2 2 -?N 
1-/3 ' 



The other cases {(3 = 1 and f3 > 1) are similar and left to the reader. Lemma [6. II is proved 
in the particular case of d = 1. 

The case of higher dimension is nothing but a technical complication of the previous 
case d = 1. Let us fix d > 1, a = (ax, ... , ad) £ 1>k and denote by j = (ji, ■ ■ ■ ,jd) any 
element of Z d . 

Let us begin with the case where a £ T>k- Let (Ji, . . . , J^) the d integers between — N 
and N such that for all I = 1, . . . , d, J\ ^ 2aiN ^ J + 1, with at least one inequality that 
is strict. The coordinates J/ and J\ + 1 are by construction the closest integers to 2a/ N 
in — N, . . . , N. For the rest of this proof, we will refer to them as critical coordinates. 
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Then, one can decompose the sum V ■ 



2N 



-p 



according to the number p of critical 



coordinates among (ji, . . . , j^) = j, where jsa typical index: 



y — -a =y t y — 

Aj 2N Aj Aj Aj 9.m 



P=0 (ii,...,ip) jEJ(i lt ...,i p ) 



2N 



(A.l) 



where the second sum is taken over all the vectors (i\,...,i p ) with strictly increasing 
indices taken among 1, . . . , d and where ■ ■ ■ , i p ) is a notation for the set of vectors 

j = (ji, . . . ,jd) such that ji x is critical for every I = 1, . . . ,p. 

In the sum (jA.ip , let us treat the cases p = and p > separately. Let us first focus 
on the case p = 0: it corresponds to vectors j without critical coordinates, which means 
that we restrict ourselves to j such that for every k = 1, . . . , d, either jf. < (in such case 
\j k - 2a k N\ = 2a k N-j k ) or either j k > J k + 1 (in such case \ j k - 2a k N\ = j k - 2a k N). In 

particular, this sum can be divided into 2 d sums XljeD 

subdomain of [—1/2, l/2] d , which is defined by this binary choice for each j k . For simplicity, 
we only treat the case of Dq := {j = (ji, . . . ,jd)\ Vfe = 1, . . . , rf, j k < J k }. The case of the 
other 2 d — 1 subdomains can be treated in a similar way. 
We have successively, 



j_ 

2N 



-P 



where D is a connected 



y — 

A^ 2N 



jeD 



2 P N P ^ 

j k <Jk-l 



3l) 



i=i 



-13/2 



<2V 



J-N i- 



Jd 
N 



r N + 2a 1 N r 
J2aiN-Ji J2 



J](2a { JV-t,) : 



-/3/2 



dti . . . dt„ 



r 2Ar i 

^ CN 13 ~^r d - 1 dr, 



-/3/2 



dui . . . dn^ 



(A.2) 
(A.3) 

(A.4) 
(A.5) 



where wn > is the distance to of the point of coordinates {2a\N — Ji, . . . , 2a^ — J^). 
The estimates found in Lemma 16.11 are then straightforward: for example in the case 
P < d, an upper bound for the last quantity is CN" N d ~P = CN d . The other cases are 
treated in the same manner and lead to the same desired estimate. 



As far as the case < p ^ d is concerned, the particular case p = d is a bit special: it 
corresponds to vectors j with only critical coordinates. Since in that case, each 



is either equal to 



2N ak 



or 



2N 



ak 



and is anyway larger than pn,k ^ 



l 



2NK 



2N ak 

(where 



the quantity pn,k has been defined in the beginning of this proof), the contribution of 



this case to the whole sum can be bounded by 2° 



^AnPrp = cnPrp. 
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Let us now concentrate on the case < p < d: Then for a fixed choice of indices 
(ii, . . . , i p ), we have 



E 

iej^,...,^) 



3 

2N 



- 2 

jGj (il ,...,i p ) 

je^i!,...,^) 



i=ii 



ll—J 

s 



f A. 

V2A^ 



l=£il,...,lp 



3i 

2N 



ij^ii,...,i p 
-P/2 



\2N 



-13/2 



-13 



But this last sum is nothing else than 2N ~ a > where a (resp. j) is the vector in 

built upon the vector a (resp. j) with all its coordinates of index in {ii, . . . , ip} 
removed. Since p > 0, we see that, by induction hypothesis, that the previous sum can be 
bounded by 

jCN d - p K d - p In N if0<d-p 
[CN? jf/3>d-p. 

In particular, if (3 ^ d, then the contribution to (jA.ip of the sum over < p < d can be 
bounded by CN d -PR d -P In N ^ min (CK d N d In N, CN?) . If /3 < d, it is also straightfor- 
ward to see that this contribution is also smaller than CN d K d . The proof of Lemma 16. II 
follows, by induction. □ 
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